SISSA 74/96/EP 



The Bosonization of the 
Electroweak Penguin Operators 



M. Fabbrichesi and E.I. Lashinjj] 

INFN, Sezione di Trieste 
and 

Scuola Internazionale Superiore di Studi Avanzati 
via Beirut 4, 1-34013 Trieste, Italy. 



Abstract 

We give the complete 0(p 2 ) bosonization of the electroweak Penguin operators Qj y s 
and compare the result with that of the gluon Penguin operators Qs^- We find 
that, in addition to the usual (constant and current-current) parts, there are three 
new terms not discussed previously in the literature. Two of these are present in the 
factorization approximation and should be included in the standard definition of the 
-E^g-factors. The impact of these corrections on the direct CP- violating parameter 
e'/e is briefly discussed. 



SISSA 74/96/EP 
May 1996 

^Permanent address: Ain Shams University, Faculty of Science, Dept. of Physics, Cairo, Egypt. 



1. Penguin operators play an important role in the physics of the standard model 
at low energies and, in particular, in the non-leptonic decays of kaons. The estimate 
of their matrix elements has been a controversial subject from the very beginning 
because of a subtle cancellation of the leading contribution in chiral perturbation 
theory. 

Most of the debate [|I|, 0, however, was centered around the gluon Penguin 
operators — at the time the only relevant ones — that we write as 

Qb = Nv-AEfe)v+A (!) 

Qa = (s a ^) v _ A ^(g^ a ) v +A (2) 
<? 

by means of the by- now-standard notation in which a, (3 denote color indices (a, (3 = 
1, . . . ,N C ), the subscripts (V ± A) refer to the chiral projections 7^(1 ± 75) and 
color indices for the color singlet operators are omitted. The other gluon Penguin 
operators ((^3,4) discussed in the literature are readily bosonized as the product of 
two left-handed currents and we will not discuss them. 
The electroweak Penguin operators f|] are defined as 
3 

Q7 = ^(^V-AE^feW ( 3 ) 
z q 
3 

Qs = ^(s«^) V -aE^(?^«)v+a, (4) 
A q 

where e q are the quark charges (e d = e s = — 1/3 and e u = 2/3). The appereance of 
these operators has not revived the discussion on their bosonization and they have 
been treated along the same lines as the gluon ones. Yet, new features emerge when 
considering this new class and the argument followed in the bosonization and the 
determination of the chiral coefficient of the gluon Penguin fail us here. Let us see 
why. 

2. In order to proceed we first write the two electroweak operators as 

3 3 
Q7 = -^e d Q 5 + -(e u - e d )AQ 7 

3 3 

Qs = ^e d Q 6 + ~(e u - e d )AQ 8 (5) 
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where 



AQ 7 = (sd) Y _ A (uu) v+A 

AQ 8 = (s a d^) y _ A (M /3 Ma) v+ A • ( 6 ) 

This splitting has the advantage of separating out the pure octet part of the elec- 
troweak operators — that is, the part already present in the case of the gluon Pen- 
guins and which therefore does not introduce any new features. 

To obtain the chiral representation of the operator AQ 7 8 we must consider the 
bosonization to 0{p 2 ) of the quark densities: 

7 2 

4 

7 2 



QxQr 



-2B, 



-E + L 5 ELLENS + AB L S ZM ] E 



<?R<ll 



~2Bn 



-E* + L 5 E t D u ED M E t + AB Q L^M E f 



(7) 



i.i 



where B = — (qq) / f 2 , M. = diag[m n , m^, m s ] and the subscript ij is the flavor 
projection. The bosonization in eq. (]?]) can be obtained by either considering directly 
the generalized mass term in the strong lagrangian or by computing 

S £XPT SC xPT 



and 



(8) 



5M^ 5M 
respectively, where £ xPT is the strong lagrangian to 0(p A ). 

The possible terms containing second derivatives have been eliminated in eq. (0) 
by means of the equations of motion. If we do not use the equations of motion and 
retain the second derivative terms there are two additional terms that in principle 
must be added 

2B a £> 2 E + 2B c 2 EZ^EE 1 . (9) 

However, only two out of the three derivative terms in eq. (0) and eq. 
independent because we can eliminate one of them by means of the relation 

ED 2 EE t + 2L>„E J D /i E t E 



arc 







(10) 



which follows by the unitarity of the matrix E alone. The addition of these terms 
does not modify our result based on eq. (|7|) except for a non-factorizable term || 



c nf Tr (A^E A^E f 



:n) 
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which is generated by the second order derivative acting on the two densities simul- 
taneously. Such a term goes beyond the vacuum saturation approximation (VSA) 
and we shall neglect it. 

Additional terms proportional to the counterterms L 4 , L 6 and L 7 are not neces- 
sary because they give contributions to eq. (0) of the form, for instance for L4: 

2B a UTt{D^D i X)T lij: (12) 

which only represent a wave-function renormalization induced by the strong sector. 

We prefer to use the equations of motion because we only need the higher order 
terms at the tree (on-shell) level and — as stressed in || — the counterterms neces- 
sary in the renormalization procedure are only correct if a classical background is 
assumed. This is the same procedure followed in the definition of the 0(p 4 ) chiral 
lagrangian C xPT on which eq. (0) is based. 

The crucial point is that — no matter what procedure one follows — there are still 
two independent terms in eq. (^) beside the leading (constant) one. 

We now turn to the bosonization of the electroweak operator Q 8 . By applying 
eq. (0) to the operator 

-l2{s L u R ){u R d L ) (13) 

which is obtained by a Fierz transformation from the operator 3AQ$/2, we obtain 
three terms 

-3((gg» 2 Tr (\^)Tt (A*e) (14) 



,{qq) 2 L 



-12"" ^ 



Tr Tr (\ l 2 D^D^ S f 



P 

+ Tr (XID^D^ E) Tr (\\T?)] (15) 



(qq) 2 BoLi 



48 N " 7 „7"" 8 [Tr (aJjS 1 ") Tr (A^EA^E 



P 

+ Tr (A?E) Tr (A^E^E^l , (16) 



where the projection matrices are defined by (Xftik = bubjk- Another term 

- ^ Tr (A^E^E) Tr (AiED^) (17) 
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is obtained by considering the leading order bosonization of the quark currents in 
the operator Q 8 as 

«£tVl - -£&D^) (18) 

- -4( e ^ e %- (i9) 

Such a term is completely determined, being of the leading order, and requires no 
further discussion. 

The current literature on the VSA estimate of the electroweak operators |4]] has 
dealt only with the terms (|14|) and fll7D . Ref. || discusses the complete lagrangian 
by means of a different bosonization technique which also include non-factorized 
configurations. In comparing our lagrangians with that written in ||, care should 
be taken in rewriting single traces as double traces and vice versa (see the appendix 
of ref. §). 

The usual expression for the bosonization of the pure octet part (the only part 
for the gluon Penguins Qs^) is readily obtained from eqs. (|T^)-([T6|) by replacing the 
projection on the u (1) quark by the sum over all quark flavors to obtain: 

P 



_ 24 {m) L " Tr (XID^D^) (20) 



-48 Tr 



\1(M*E + T!m)] , (21) 



where the leading constant term has vanished by projection and the mass correction 
is just a renormalization that can be reabsorbed in the definition of x = ^B^M. in 
the strong lagrangian. 

In this case there is only one possible bosonization term ( p0| ) and one constant 
(L 5 ) to be determined. The usual bosonization for Q 6 (and eq. (^) below) is 
obtained by multiplying eq. (EDI) by 2/3. 



3. If we neglect the non-factorized term([TT|) and the mass term (16, then the 
chiral lagrangian necessary in the VSA approximation is determined by the same 
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argument used in the case of the gluon Penguin. In particular, it is possible to fix 
L 5 by computing the ratio between the kaon and pion decay constants: 

A(K + -> (i+uj 



which gives, in the large 1/N C limit [0, 



(22) 



4 V U J m K - m l 



L 5 = ~ ( 1 - ¥ I =5 = 2 (23) 



and 

L 5 (m p ) = (1.4 ±0.5) 10" 3 , (24) 

if chiral logarithms are kept 0. 

In order to include the mass-term correction, the constant L 8 can be determined 
from the GMO formula and the knowledge of L 5 . It is found H that 

L 8 (m p ) = (0.9 ±0.3) 10~ 3 . (25) 

We now write the matrix elements of the operator Q 6 and Qg (those of and 
Qi are readily obtained by going to the next order in 1/N C ) for the decay of a kaon 
into two pions. As usual, we split them into isospin amplitudes (I = 0, 2) and find 

(Q«f = ( 26 ) 



where X = \f%f % {rr? K — ml). Eq. (p6|) is the usual VSA result for the gluon Penguin 
operator. 

The operator Q 8 dominates the channel 1 = 2 where we find that 
(Qs)^ = v/6 ± 4v/6 (4L 8 - L 5 ) m\ - ^-X (27) 



The second term in eq. ( [27D is the result of the two new terms we are discussing. It 
is as large as 20% of the leading one, as opposed to the small momentum correction 
contained in the third and last term that is only 1%. 
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As we have already pointed out, in the literature 0, only the leading (con- 
stant) term and the X-term have been so far included. As a consequence, the usual 
definition of the B$-f actor 

I r~\ \ an y model 

\ Qb) 

(QsY 

that quantifies deviation from the VSA is not correct and should be changed ac- 
cording to eq. Q27D- 



B * = TTTTvsa ' ( 2J 



5. Because of the importance of electroweak operators in the determination of 
the direct CP- violating parameter e'/e, we have estimated their effect in the VSA 
and compared the result that includes the new terms with that that does not. As it 
can be seen in Fig. 1, the effect is about 20% in the range of input parameters we 
considered. We should however bear in mind that the cancellation between gluon 

878 
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-0.001 



Q3 Q4 Q5 Q6 Q7 Q8 Q9 Q10 All 

Figure 1: Contributions of the various operators in the VSA to e'/e at fi = 0.8 GeV. 
Grey (black) bars (do not) include the extra terms. All operators are shown for 
completeness. 

and electroweak operators is not as effective in the VSA as it is in all more refined 
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computations f?|, [|, as shown by the rather large final value of e'/e in Fig.l. For 
this reason, in any computation in which such a cancellation is more complete, the 
impact of the new terms could be much more dramatic 



6. For the sake of comparison, we now turn to a specific example of modeling of 
low-energy QCD. In the chiral quark model (xQM) (see ref. f| for discussion and a 
list of references) we have 



-xQM 



txQM 
^8 



8M| (qq) | 
16M\ (qq) 
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Ml) 

Mf_ 

I (qq) 



1 N r 



24 16tt 2 



where M is a parameter characteristic of the model. 
We thus find that 



(Qs) 



xQM 

2 



V6 



(qq)' 



(qq? 



c n{ m 



2 + ( 4L xQM 



f3 v - f5 

where also the non-factorizable term proportional to 



L xQM\ 



m 



K 



1 

C n f = ~ 



4 M\ (qq) \ 

is included for completeness. 

In particular, the combination 

3 / 4 



4L 



xQM 



-xQM 



3 / 4 M 



M| (qq) | 4 A 2 1 (qq) 



1 

4 



1 N r 



(qq) 



8 16tt 2 



(29) 
(30) 



'2N r 



X,(31) 



(32) 



(33) 



turns out to be numerically of the same order as the VSA result. An estimate of 
e' I e in this model (albeit in the chiral limit inclusive of only the second of the three 
terms in eq. (EJ3I) and of the one proportional to c n f) is presented in ref. M. 



M.F. would like to thank E. de Rafael for discussions. 
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